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Abstract 

A two-stage model reduction methodology, combining the classical Component Mode 
Synthesis (CMS) method and the newly developed Enhanced Projection and Assembly 
(EP&A) method, is proposed in this research. The first stage of this methodology, called 
the COmponent Modes Projection and Assembly model REduction (COMPARE) method, 
involves the generation of CMS mode sets, such as the MacNeal-Rubin mode sets. These 
mode sets are then used to reduce the order of each component model in the Rayleigh-Ritz 
sense. The resultant component models are then combined to generate reduced-order system 
models at various system configurations. A composite mode set which retains important 
system modes at all system configurations is then selected from these reduced-order system 
models. In the second stage, the EP&A model reduction method is employed to reduce 
further the order of the system model generated in the first stage. The effectiveness of the 
COMPARE methodology has been successfully demonstrated on a high-order, finite-element 
model of the cruise-configured Galileo spacecraft. 


1. Background and Motivation 

Multibody dynamics simulation packages are gaining in popularity among dynamists 
for the simulation and analysis of systems of interconnected bodies (some or all of which are 
flexible). One such program, DISCOS , 1 was used in the development of control systems on 
board the Galileo spacecraft. The dual-spin Galileo spacecraft was modeled as a three-body 

55 

PRECEDING PAGE BLANK NOT FUMED 


system, consisting of a flexible spinning rotor, a flexible stator, and a rigid scan platform. 

For complex systems such as the Galileo spacecraft, practical considerations (e.g., sim- 
ulation time) impose limits on the number of modes that each flexible body can retain in a 
given simulation. Modal truncation procedures must be used to select and retain a limited 
number of “important” modes which capture the salient features of the component dynam- 
ics. The Enhanced Projection and Assembly (EP&A) 6 ’ 7 technique is one way of performing 
this task. 

The EP&A method, 6,7 is a model reduction methodology for articulated, multi-flexible 
body systems. In this method, a composite mode set, consisting of “important” system 
modes from all system configurations of interest, and not just from one particular system 
configuration, is first selected. It is then augmented with static correction modes before be- 
ing “projected” onto the component models to generate reduced-order component models. 
To generate the composite mode set, eigenvalue problems concerning the full-order system 
models, at all configurations of interest, must be solved repetitively. This is a drawback of 
the EP&A method because solving large eigenvalue problems can be costly. To overcome 
this difficulty, a two-stage model reduction methodology, combining the classical Compo- 
nent Mode Synthesis (CMS) method and the Enhanced Projection and Assembly method 
(EP&A), is proposed in this research. 

The stages involved in the proposed technique, to be called the COmponent Modes 
Projection and Assembly model REduction (COMPARE) method, are illustrated in Fig. 
1. First, CMS mode sets, such as the MacNeal- Rubin mode sets, are generated and used 
to reduce the order of each component model in the Rayleigh-Ritz sense. These compo- 
nent mode sets are then assembled using the interface compatibility conditions to generate 
reduced-order system models at various system configurations. The order of these reduced- 
order system models is typically smaller than that of the full-order system model. 

In the second stage, the newly developed EP&iA model reduction method is employed to 
reduce further the order of the system model generated in the first stage. As described above, 
the composite mode set is augmented with static correction modes before being projected 
on the CMS-generated component models to generate the final reduced-order system model. 
In this way, COMPARE retains the merits of both the CMS and EP&A methods, without 
their demerits. The effectiveness of COMPARE will be verified using a cruise- configured 
Galileo spacecraft model. 
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2. Component Mode Synthesis Method Revisited 2,8 

The Component Mode Synthesis (CMS) method is a Rayleigh-Ritz based approximation 
method that is commonly used to analyze linear, high-order structural dynamics problems. 
To use this method, the structure is first subdivided into a number of components (or 
substructures), and a Ritz transformation is employed to reduce the model orders of these 
substructures. Many component mode sets may be used to perform this reduction but 
the MacNeal-Rubin (M-R) and Craig-Bampton (C-B) mode sets were shown to have good 
convergence properties in the sense of CMS. Once reduced, the reduced-order component 
models are then coupled using the interface compatibility conditions to form the reduced- 
order system model. 

Since the Craig-Bampton and MacNeal-Rubin mode sets will be used in the present 
research to reduce the orders of the component models, their constructions are first briefly 
reviewed. To this end, consider a multi-flexible body structure as depicted in Fig. 2. 
The undamped motion of each component of the structure can be described by a matrix 
differential equation 

Mnnin “b Ffi 1 

where is an n x 1 displacement vector, and M nn and K nn are n x n mass and stiffness 
matrices of the component, respectively. Note that the matrix dimensions are indicated by 
the matrix subscripts. The n x 1 force vector acting on the component is denoted by F n . A 
similar equation can also be written for component B. 

To generate the C-B or M-R mode set, the last equation is partitioned as follows . 

FI ij 
Mji Mjj 

where x, and xj represent the interface and interior coordinates, respectively. 

2.1 Craig-Bampton Mode Set 

The Craig-Bampton mode set is generated by augmenting a low-frequency subset of the 
fixed interface (I/F) normal modes with a set of static-shape functions termed constraint 
modes. The first k fixed I/F normal modes $,■*, and the ordered eigenvalue matrix Ajtjt are 
related by the following relation : 

—Mjj$jkhkk + Kjj$jk = 0 - ( 3 ) 
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There are ways to decide on the number of modes to be kept in One way is to 

keep all modes whose frequencies are less than twice a characteristic frequency of the system 
(e.g., control bandwidth). 2 The above determined normal modes are then augmented with i 
constraint modes, where constraint modes are static-shape functions that result by imposing 
unit displacement on one coordinate of the ?-set while holding the remaining coordinate in 
that set fixed. It can be shown that the interior displacement for these constraint modes is 
given by 8 

*ji = -A'" 1 Kn . (4) 

The C-B mode set is 

Oik lit 

which is then used to reduce the full-order component model. 

2.2 MacNeal-Rubin Mode Set 

In a manner similar to generating the Craig-Bampton mode set, the MacNeal-Rubin 
mode set is generated by augmenting a low-frequency subset of the free I/F normal modes 
with a set of static force response functions termed residual modes. The first k free I/F 
normal modes and the ordered eigenvalue matrix A **, are related by the following 

relation : 

A/ nn $ n Jt Afcfc "4* A nn $ n Jt = 0 • (b) 

The kept eigenvector matrix which has been normalized with respect to the mass 
matrix, may be partitioned into its rigid-body and flexible parts: [<£ nr <!>„/]. Let A// be the 
eigenvalue matrix associated with the kept flexible modes. Then, the residual modes may 
be determined by 

- (Pnn Snn Pnn ~ *$n/ A fj ) Pna > (A) 

where P nn = — M nn <£ nr <$J r , 5„„ is the “pseudo” flexibility matrix of the component, 8 

defined as follows. Let the set of all physical coordinates be divided into three subsets : r, a, 
and w. The r set may be any statically determinate constraint set which provides restraint 
against rigid-body motion. The a set consists of coordinates where unit forces are to be 
applied to define attachment modes (i.e., at the interface coordinates, and at coordinates 
where external forces are applied). Finally, the w set consists of the remaining coordinates 
in z. Using these definitions, the component stiffness matrix K nn is partitioned as follows : 
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The pseudo-flexibility matrix 5„„ is now given by a matrix of the form 8 
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Finally, the matrix F na is given by [0 ar , Iaa , 0 au ,] T , where the identity matrix I aa is asso- 
ciated with the a interface coordinates. The M-R mode set is [$„* 'Fna], which is then used 
to generate a reduced-order model for the component. 

3. A Component Modes Projection and Assembly Model Reduction (COM- 
PARE) Methodology 

Once CMS-based reduced-order component models are generated, they are assembled 
using the interface compatibility conditions to produce reduced-order system models at vari- 
ous system configurations of interest. Since the orders of these reduced-order system models 
are typically smaller than those of the full-order system models, we have accomplished the 
first of the two model reduction steps of the COMPARE methodology. However, note that 
these CMS-generated reduced-order system models were obtained without using knowledge 
of any system-level input-output information. This drawback is remedied in the second 
stage, in which the EP&A methodology 6,7 is used to further reduce the order of the models 
generated in the first stage. Since the EP&A methodology has been described elsewhere, 6 ’ 7 
it will only be briefly reviewed here. 

Consider a system with two flexible components. The undamped motion of component 
A, as described by either its M-R or C-B mode set, is given by 

Ippiip + AppT?p = G pa u a , 

vt = • 

Here, t) p and A pp are the generalized coordinates and the diagonal stiffness matrix of com- 
ponent A, respectively. The dimension of rj p is p. The matrix G pa is an control distribution 
matrix, and u* is an a x 1 control vector. Similarly, the matrix H* p is an output distri- 
bution matrix, and is an b x 1 output vector. Similar equations can also be written for 
component B. 
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The system equations of motion at a particular articulation angle a may be constructed 
using these component equations, and enforcing displacement compatibility conditions at 
the component interface. To this end, let P(a) =[P^ (<^), Pq C (°0] T he any full-rank matrix 
mapping a minimal system state rj e into 




'P p A MY 

U?J 


e®(a). 


[ J ?e] » 


( 8 ) 


where r) e is an e x 1 reduced-order system coordinate, and e = p + q — i {i is the number of 
I/F constraint relations). For ease of notation, the dependencies of the matrices P A e , P®, 
r/p , y®, and t] e on a are dropped in the sequel. Substituting rj A = P A t rj e and r/® = P® r) e 
into (7) and the corresponding equations for component B, pre- multiplying the resultant 
equations by P A * and P® > respectively, and summing the resultant equations gives 


Afee^/t 4" Pce'He G^a^a , 

Us H 3e Tj e , 


(8a) 

m 


where M ee , K ee , G ea , and H ae , all functions of a, are given by 

M — P aT P A 4- pB T p B 
K - P aT A a P a 4- P® T A®P® 

iV ee — 1 pe il pp J pc * 2 qe lv qq 2 qe ’ 

n _ pA T A . pB T (~<B 
— 1 pe ^ pa ' 1 qe ^qa ’ 

where y a = [ y AT y® T ] r , and s — b + /. To arrive at the equation for G ea , we have as- 
sumed that u A — u ® = u a . Otherwise, the term G ea u a in (8a) should be replaced by 
\pA qA pB qB i r u a U B i T gj nce \ s symmetric and positive-definite while I\ ee is 

symmetric and positive semi-definite, a transformation $ ce that diagonalizes M ee and I\. ee 
simultaneously can always be found. Let <f> e be the corresponding generalized coordinate, 
i.e., r/e = $ee <t>e • Substituting this relation into (8a), and pre- multiplying the resultant 
equation by gives 


r i A r>A 
12 bp 2 pe 

i jB jdB 

" la 2 at J 


fee ^tt$t — 5 . . 

v”) 

Vs = H st <f>e 1 

where $J e M ee $ee — he, Hse — Hse$ee, G ea = $J e G ea , and where Aee — $J e hee$ee 
contains the undamped, reduced-order system eigenvalues along its diagonal. Equation (9) 
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represents the reduced-order system model obtained from the first stage of the COMPARE 
methodology. 

The EP&A method is used in the second stage of COMPARE. With the EP&A method, 
only k of the system’s e modes are kept while the remaining t (= e - k) modes are removed. 
The kept mode set is a composite mode set, consisting of “important” system modes from 
all system configurations of interest, and not just from one particular configuration. With 
this understanding, we have 


77e = $ee<£e = [ $ek 



= $ek4>k , 


( 10 ) 


where <f> k and <f>t are generalized coordinates associated with the “kept’ and truncated 
modes, respectively, and and $ e t the corresponding eigenvector matrices. 


The composite mode set $ e Jt may now be projected onto the CMS-generated compo- 
nent models: = P £ $ ek <t> k = ■ Here > $ k denotes reduced sets of generalized 

coordinates of component A. The substitution of the last relation into (7) produces the 
“constrained” equations of motion for component A: 


< *ik it + < a t P *tk rt = << g 
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( 11 ) 


A second reduced-order system model can now be constructed using (11), a similar 
equation for component B, and the displacement compatibility conditions at the component 
I/F. The order of this new reduced-order model is smaller than that obtained from the first 
step (cf. (9)) due to the truncation of “t” modes in (10). In addition, it has been proven 
that the modes retained in the composite mode set $ efc are captured exactly by the resultant 
reduced-order system model (with a number of extraneous modes > > ). However, the 
static gain of the resultant reduced-order system model is not the same as that given by 

(9)- 7 


Two different approaches were introduced in Lee and Tsuha 7 to preserve the static 
gain of the system described by (9) in the second reduced-order system model. The first 
approach involves augmenting the k “kept” modes of the system with an additional a modes 
so as to create a statically complete mode set. The enlarged mode set is then “projected” 
onto the components, and the reduced components are assembled as usual. 7 In the second 
approach, the k “kept” modes of the system are first projected onto the components. After 
the projections, the reduced-component mode sets are each augmented with static correction 
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modes. The augmented mode sets are then used to reduce the components, and the resultant 
component models are assembled to generate the reduced-order system model. The details 
of these approaches are given below. 


3.1 Component-level Augmentation Techniques 6,7 


In the component-level augmentation approach, the k “kept” modes are first projected 
onto the components. We first write 


Ppr r*,\ vi 


k > 


Vp — 'bpk —kk u k — ~fpk v k 
1/f = «?* S& uf = Tf t u? = [Tf r T?,] u* , 


( 12 ) 


where E** is the eigenvector matrix associated with an eigenvalue problem of (11), and u k 
is the corresponding generalized coordinate. The transformation matrices T pk and are 
defined in (12). Matrix partitions T^ r and T p j denote eigenvectors associated with the 
rigid-body and flexible modes of the projected component A model. Similar partitions are 
made for the eigenvectors associated with the projected component B model. 

Let Tp r , T etc. be normalized such that 

T f T Tp r = I rr , T f r App Yp r = 0 rr , 

rff = Iff , T f } App Tpy = A?, , 


(13) 


where A^j is an eigenvalue matrix associated with the projected flexible modes of component 
A. Similar expressions can also be written for component B. Using the matrices defined in 
(12-13), residual modes, 8 described in Section 2.2, may once again be used to augment the 
projected mode sets of the components. From (6), the residual modes are given by 


T * = (P£ SI P. 
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(14) 

'yA *rA T 
A pr pr ? 

(15) 


and S pp is the “pseudo” flexibility matrix of component A. If we assume that the diagonal 


stiffness matrix of component A (see A pp in (7)) is ordered so that all the rigid-body modes 
(with zero frequency) are given first, then S pp is given by 


0 0 
0 App 1 


( 16 ) 
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Here p = p - r is the total number of flexible modes in the CMS mode set of component A. 
In (14), the matrix F ' A is given by [0 ar , J««, 0 ap -r-«] T . Similar expressions can also be 
written for component B. 


The projected mode sets T^ k and T^. 
modes 


if = [Tf* TfJ 

if = prft rfj 


can now be supplemented with the residual 


f 


= T- 4 i/ A , 

— ± qv 5 


(17) 


where v A and uf are generalized coordinates associated with the residual modes. Using the 
Ritz transformations suggested in (17), the new component-projected equations of motion 
are 
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(18) 

(19) 


Using (18-19), the reduced-order system equations of motion at a particular system 
configuration a can now be formed by enforcing displacement compatibility at the interface 


C.»lf + C?(a)lf = 0 , (20) 

where C, A and C\ D q are matrices that establish the constraint relations between the gener- 
alized coordinates of CMS generated component models, and i is the number of constraint 
relations. Using (17), (20) becomes 


f r*A*rA s-*B~fB ] 
V-'ip 1 pv ^ iq 1 qv I 


v 

,B 


= [D ie ] 


= 0, 


( 21 ) 


where [D ic ] is defined in (21), and c = 2v. To construct the reduced-order system model, we 
partition the “compatibility” matrix [D IC ] using the Singular Value Decomposition (SVD) 
technique 


v;T 

v 2 


[Ac] = \C* T£ , Cf q T® ] , = [Uii} [E« , O ld ) 

where d = c - i, and [ O ld ] is an i x d null matrix. We can now write 3 ' 6,7 


( 22 ) 
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where [-/^crf] —- [V c d] (cf* (22)) is a full column rank mapping matrix, and v d denotes a 
minimum set of generalized coordinates of a statically complete reduced-order system. Sub- 
stituting v = [P A d ] v d and v* = [P? d ]v d into (18) and (19), pre-multiplying the resultant 
equations by P A J and , respectively, and summing the resultant equations give 


M dd v d + Kddv d = G da u a , 
2/s = H sd v d , 


(24) 

(25) 


where M dd , K dd , G da , and H sd , all functions of a, are given by 
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(26) 

(27) 

(28) 

(29) 


Since Mjd is symmetric and positive-definite while K dd is symmetric and positive semi- 
definite, a transformation 4?dd that diagonalizes M dd and I\dd simultaneously can always 
be found. Let $ d d be normalized with respect to the mass matrix, and let Xd be the 
corresponding generalized coordinate, i.e., 


Vd = $dd Xd ■ 


(30) 


Substituting (30) into (24-25) and pre-multiplying the resultant equation by <fr dd give 


IddXd "4" K d dXd — G da U a , 

2 Is = H sd T} d , 


(31) 

(32) 


where 


H ad = H sd $dd, (33) 

Gda = ^ d Gda, (34) 

Add = $ dd K d d$dd ■ (35) 


Here A dd is a diagonal matrix with the undamped, reduced-order system squared frequencies 
along its diagonal. It was proven in Lee and Tsuha^ that &re captured exactly in *f* d d 
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despite the augmentations of the projected component models with residual modes. The 
matrix Add also contains a number of extraneous modes. 6 


3,2 System Level Augmentation 


(36) 


The equation (9) may be decomposed into its kept and truncated parts. 

Ikk4>k + Akk<f>k = $Jk &ea U a , 

+ Aftfit = u a • 

The mode set $ et , which is truncated, will now be replaced by a smaller but statically 
equivalent mode set $ ea - To find consider the following Ritz tiansformation 

(f>t — Rta$a ? 

where R ta = A^ 1 G ea • 

Here <fi a is the generalized coordinate associated with the augmented mode set $ e <z- It 
can be shown that the static gain due to the mode set [$ c fc ^eo ] = [$efc QetRta] is identical 
to that of the original system (cf. (9)). 6 Hence, [$ e fc ^et-Rta] is a statically complete mode 
set. 


(37) 


Several observations regarding the augmented mode set [4>eJt 4> e t R«a] are in order. First, 
note that the augmented mode set contains two parts. The first part contains a selected 
number of eigenvectors which satisfy the eigenvalue problem defined by M ee and K ee . The 
second part is formed using the residual flexibility matrix and the distribution matrix of the 
external load. Hence, Rta , unlike $ e *, is a function of the external load distribution. 

Next, we observe that $ e t Rta is identical to the residual mode defined in Section 2.2. 
Also, note that this mode set is mass and stiffness orthogonal to the retained mode set $ e *. 
Hence, the two parts of the augmented mode set are linearly independent. This statically 
complete mode set is now ready for projections onto the various flexible component models 
to generate reduced-order component models. The reduced-order component models may 
then be assembled to generate the reduced-order system model using the procedures outlined 
in Section 3.1. 


4. Applying COMPARE on A High-order Finite-element Galileo Model 

The effectiveness of the proposed COMPARE methodology will now be demonstrated 
using a high-order finite- element model of the cruise-configured Galileo spacecraft. The 
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three-body topology of the dual-spin Galileo spacecraft is illustrated in Fig. 3. 7 The rotor 
is the largest and most flexible component represented, with 243 dof. The smaller and 
more rigid stator is represented with 57 dof. Lastly, the scan platform is the smallest body 
idealized as rigid, with 6 dof. 

For the purpose of controller design, a low-order system model, accurate at all config- 
urations of interest and over a frequency range of interest (0-10 Hz) is needed. To this end, 
we apply the MacNeal-Rubin version of the COMPARE methodology on the Galileo model. 
The first stage of COMPARE requires the generation of M-R mode sets for all the flexible 
components. Following standard procedures, free interface normal modes of both the rotor 
and stator are first determined, and then truncated at twice the frequency of interest (20 
Hz). Next, these truncated normal mode sets are each augmented with residual modes to 
generate the needed M-R mode sets for both the stator and rotor. 

Next, the MacNeal-Rubin mode sets and the interface compatibility conditions are used 
to construct system models at all system configurations of interest, and determine from them 
important system-level modes at all clock angles of interest. The selected composite mode 
set has 8 rigid-body and 21 flexible modes. 

The next step is to augment the composite mode set with one or more static-correction 
modes. For the Galileo example, we augment the composite mode set with two residual 
modes, one for an input torque about the Z-axis on the rotor side of the rotor/stator interface, 
and a second equal and opposite torque on the stator side of the interface. The enlarged 
mode set is then projected onto the flexible components. The resultant reduced-order models 
of the rotor and stator have 29 (with 6 rigid-body) and 21 (with 8 rigid-body) modes, 
respectively. The assembled reduced-order model has 44 (with 8 rigid-body) modes. The 
natural frequencies of these reduced-order component models and of the system model at a 
clock angle of 300 degrees are tabulated in Table 1. All system flexible modes retained in 
the composite mode set have been captured exactly in the reduced-order system model. 

Comparisons of Bode plots of full-order and reduced order models, at clock angles of 
60, ISO, and 300 degrees, are given in Figs. 4, 5, and 6 respectively. Actuation was done 
at the Spin Bearing Assembly (SBA) located at the rotor-stator interface (along the Z-axis, 
cf. Fig. 3), and sensing was done by a gyroscope located on the scan platform. From 
these comparisons, we observe that the frequency responses of the full-order models, at all 
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configurations of interest, have been closely captured by their reduced-order counterparts, 
over the frequency range of interest (0-10 Hz). Results obtained at other clock angles are 
similar to those depicted in Figs. 4-6. 

5. Concluding Remarks 

A two-stage model reduction methodology, called COMPARE, is proposed in this re- 
search. The first stage of this methodology involves the generation of CMS mode sets for 
the flexible components. The resultant component models are then combined to generate 
reduced-order system models at various system configurations. A composite mode set which 
retains important system modes at all system configurations is then determined from these 
reduced order system models. In the second stage, the EP&A model reduction method is 
employed to reduce further the order of the system model generated in the first stage. 

The merit of the COMPARE methodology is that system models (at various system 
configurations) assembled using CMS-generated component models are smaller in size than 
the full-order system models. Hence, COMPARE alleviates the need to solve large-order 
eigenvalue problems repetitively. The need to generate the components’ M-R or C-B mode 
sets is not a disadvantage because efficient software exists for their construction (see, e.g., 
Tsuha 9 ). The effectiveness of COMPARE, using the M-R version of COMPARE, has been 
successfully demonstrated on a high-order, finite-element model of the cruise- configured 
Galileo spacecraft. 

6. References 

1. Bodley, C.S., Devers, A.D., Park, A.C., and Frisch, H.P., “A Digital Computer Pro- 
gram for the Dynamic Interaction Simulation of Controls and Structure (DISCOS), 
NASA Technical Paper 1219, Vols. I and II, NASA Center for AeroSpace Information, 
Baltimore, Maryland, May, 1978. 

2. Spanos, J.T. and Tsuha, W., “Selection of Component Modes for the Simulation of 
Flexible Multibody Spacecraft,” Journal of Guidance, Navigation, and Control, Vol. 
14, No. 2, pp. 278-286, March/ April, 1991. 

3. Bernard, D., “Projection and Assembly Method for Multibody Component Model Re- 
duction,” Journal of Guidance, Control, and Dynamics, Vol. 13, No. 5, Septem- 
ber/October, 1990. 


67 



4. Eke, F.O. and Man, G.K., “Model Reduction in the Simulation of Interconnected Flex- 
ible Bodies,” paper AAS 87-455, AAS/AIAA Astrodynamics Specialist Conference, 
Kali spell, Montana, August, 1987. 

5. Tsuha, W. and Spanos, J.T., “Reduced Order Component Modes for Flexible Multibody 
Dynamics Simulations,” paper presented at the AIAA Aerospace Sciences Meeting, 
Reno, Nevada, January, 1990. 

6. Lee, A.Y. and Tsuha, W.S., “Applying the Enhanced Projection and Assembly Model 
Reduction Methodology on Articulated, Multi-flexible Body Structures,” Proceedings 
of the IEEE Singapore International Conference on Intelligent Control and Instrumen- 
tation, February 17-21, 1992. 

7. Lee, A.Y. and Tsuha, W.S., “An Enhanced Projection and Assembly Model Reduction 
Methodology,” Proceedings of the AIAA Guidance, Navigation, and Control Confer- 
ence, New Orleans, Louisiana, August 12-14, 1991. 

8. Craig, R.R., Jr., Structural Dynamics: An Introduction to Computer Methods, John 
Wiley and Sons, Inc., New York, 1981. 

9. Tsuha, W.S., The Benfield-Hruda, MacNeal-Rubin, and Craig-Bampton Mode Set Genera- 
tion Procedures, JPL Publication D-7213 (internal document), Jet Propulsion Laboratory, 
California Institute of Technology, February, 1990. 

7. Acknowledgments 

The research described in this paper was conducted at the Jet Propulsion Laboratory, 
California Institute of Technology, under a contract with the National Aeronautics and 
Space Administration. The authors wish to thank Dr. D. Eldred, Dr. J. Spanos, and Dr. 
M. Wette for many helpful discussions and valuable suggestions. We also wish to thank Dr. 
M. Lou and Dr. G. Man for their interest and encouragement. 


68 



Table 1 Frequencies of Reduced-order Stator, Rotor, 
and System (at a Clock Angle of 300 degrees) Flexible Modes 


Flexible Mode 

^stator (Hz) 

^rotor (Hz) 

^system (Hz) 

. 1 

7.105 

0.143 

0.127 

2 

9.129 

0.866 

0.864* 

3 

10.561 

1.237 

1.236* 

4 

14.560 

1.483 

1.479* 

5 

43.172 

1.728 

1.707* 

6 

50.070 

2.286 

1.734* 

7 

63.530 

2.809 

2.072* 

8 

80.867 

3.647 

2.351* 

9 

86.989 

3.996 j 

2.815* 

10 

96.605 

5.207 

3.707* 

11 

166.34 

5.337 

4.167 

12 

240.18 

5.994 

5.231* 

13 

254.52 

6.410 

5.433 

14 


9.503 H 

5.467* 

15 


10.291 

6.150 

16 


10.553 

6.436 

17 


13.536 

7.056* 

18 


15.613 

8.153* 

19 


29.188 

9.606* 

20 


41.181 

9.613 

21 


58.524 

10.294 

22 


69.003 

10.420* 

23 


77.722 

10.555* 

24 



13.534* 

25 



13.997 

26 



15.860 

27 



16.800* 

28 



20.591 

29 



21.280* 




28.462 

31 



34.316 

32 



41.214 

33 



48.011 

34 



58.318 

35 




36 



82.129' 


f Exactly captured system flexible modes. 
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o x i, Interface Coordinates 



• Xj , Interior Coordinates 


Fig. 2 Coordinate Definitions of A 2-Component System 
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Fig. 3 Galileo Spacecraft Cruise Model 2,6,7 
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Fig. 4 Bode Plot Comparison of Full- and Reduced-order Models 
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Reduced-order Model 






Fig. 5 Bode Plot Comparison of Full- and Reduced-order Models 
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Fig. 6 Bode Plot Comparison of Full- and Reduced-order Models 
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